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Abstract 

We introduce a generalized ensemble of nonhermitian matrices interpolating between the 
Gaussian Unitary Ensemble, the Ginibre ensemble and the Poisson ensemble. The joint 
eigenvalue distribution of this model is obtained by means of an extension of the Itzykson- 
Zuber formula to general complex matrices. Its correlation functions are studied both in 
the case of weak nonhermiticity and in the case of strong nonhermiticity. In the weak 
nonhermiticity limit we show that the spectral correlations in the bulk of the spectrum 
display critical statistics: the asymptotic linear behavior of the number variance is already 
approached for energy differences of the order of the eigenvalue spacing. To lowest order, 
its slope does not depend on the degree of nonhermiticity. Close the edge, the spectral 
correlations are similar to the Hermitian case. In the strong nonhermiticity limit the 
crossover behavior from the Ginibre ensemble to the Poisson ensemble first appears close 
to the surface of the spectrum. Our model may be relevant for the description of the 
spectral correlations of an open disordered system close to an Anderson transition. 
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1 Introduction 



Nonhermitian Random Matrix Models were first introduced by Ginibre in 1965 His 
motivation was to describe the statistical properties of nuclear resonances with a finite 
width in complete analogy with the description of the position of resonances by means 
Hermitian Random Matrix Ensembles as introduced by Wigner and Dyson |2j. Since 
then, eigenvalues of nonhermitian operators occurring in many different fields have been 
analyzed in terms of nonhermitian random matrix models, usually with additional ingre- 
dients. We mention several examples. The statistical properties of the poles of ^-matrices 
have been analyzed in great detail in || f|, f|. In QCD, the Euclidean Dirac operator 
in QCD at nonzero chemical potential (which can be interpreted as an imaginary vector 
potential), is nonhermitian resulting in the failure of the quenched approximation ||. 
Both this failure and the generic properties of the complex Dirac spectrum have been 
explained fully in terms of a nonhermitian Random Matrix Model with the global sym- 
metries of QCD |7|, H [l(| |, [11]]. Recently, a derealization transition was found in a 



one-dimensional lattice model with an imaginary vector potential [|T^, [13| . Statistical cor- 
relations predicted by the Ginibre ensemble have been found dissipative quantum maps 



14] , jig , fLBjl . Eigenvalue spacings of the Floquet matrix of a Fokker Planck equation have 
been described in terms of Ginibre statistics [fT7| . In [|Tj| an ensemble of asymmetric 
real matrices, closely related to the Ginibre ensemble, was utilized to model the dynamics 
of a neural network. 

Among more mathematically oriented works we mention the exact calculation of the 
correlation functions of an ensemble of normal random matrices with an arbitrary polyno- 
mial probability potential pO , 21 1 . Nonhermitian ensembles have been analyzed in terms 
of associated hermitian ensembles [22, [25| . Correlations of eigenfunctions have been stud- 
ied in the Ginibre ensemble [E4| . Another intriguing application is the description of an 
analytic curve by the boundary of the support of the complex spectrum of a nonhermitian 
Random Matrix Theory [^, |26|. Finally, we point out that there are interesting relations 
between the eigenvalues of complex matrices and the positions of particles in certain two 
dimensional physical systems ]2S|, [27], P9"| . For example, the Ginibre model is equivalent 
to a Coulomb problem in two dimensions [l|]. 

Based on the magnitude of the imaginary part of the eigenvalues we distinguish two 
types of nonhermiticity: weak nonhermiticity and strong nonhermiticity. Weak nonher- 
miticity is the limit of large matrices when the imaginary part of the eigenvalues remains 
comparable with the mean separation of eigenvalues along the real axis. This limit was 
identified in [|30], |3T|, [J2j, but was used earlier in the statistical theory of ^-matrices ||. 
Strong nonhermiticity refers to cases for which the real and imaginary parts of the eigen- 
values remain of the same order of magnitude in the thermodynamic limit. In this article 
we consider both types of nonhermiticity. 
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An important concept in the understanding of disordered systems is the Thouless en- 
ergy. We will define this energy scale as the energy difference below which the eigenvalues 
are correlated according to Random Matrix Theory. In diffusive disordered systems, in 
the thermodynamic limit, both the eigenvalue spacing and the Thouless energy approach 
zero whereas the number of eigenvalues in between them approaches infinity. In this arti- 



cle we will consider critical statistics |33| , |34| . |35| . |36| which refers to the case that the ratio 



of the Thouless energy and the eigenvalue spacing remains finite in the thermodynamic 
limit. A Hermitian Random Matrix model for critical statistics was proposed in |3?| . 



In that model the correlations of the eigenvalues decay exponentially beyond a Thouless 
energy resulting in an asymptotically linear behavior of the number variance with slope 
(level compressibility) less than one. In this article we generalize this model to complex 
eigenvalues and analyze its properties. In the Ginibre model the two-point correlation 
function of eigenvalues in the bulk of the spectrum drops off exponentially on the scale 
of the distance between the eigenvalues. It is therefore no surprise that we will find the 
same bulk correlations in such generalized Ginibre model. However, we find nontrivial 
long range surface correlations, characteristic of a two-dimensional Coulomb liquid. In the 
case of weak nonhermiticity we expect to find critical statistics similar to the Hermitian 
model. The analysis of this case is the main objective of this article. 

Critical statistics is associated with the multifractal behavior of the eigenfunctions 



36, p^ , |39| . The critical Hermitian model introduced in [[37]] has the unitary invariance of 
the Gaussian Unitary Ensemble with eigenvectors that are distributed according to the 
measure of the unitary group. This is no contradiction: multifractality of wave functions 
occurs in a specific basis in which disorder competes with a hopping term. Indeed, in 
BO, [4"T| it was found that the fractal dimension of the wave function determines the 



asymptotic slope of the number variance. 

Among others, critical statistics have been utilized to describe the spectral correla- 
tions of disordered system at the Anderson transition in three dimensions []3~3| , |42|| , two 
dimensional Dirac fermions in a random potential H3|, quantum Hall transition MM and 



QCD Dirac operator in a liquid of instantons Ej| |4q| . The scope of universality of critical 
statistics is still under debate. 

Our Random Matrix Model is introduced in section 2. The case or strong nonher- 
miticity and weak nonhermiticity are analyzed in sections 3 and 4, respectively. Among 
others we derive a closed expression for the two-point correlation function in both limits. 
Results for the number variance are discussed in section 5 and concluding remarks are 
macritical statistics is still under debate.de in section 6. 
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2 Introduction of the model 



Recently, a Hermitian random matrix model for critical statistics was introduced by 
Moshe, Neuberger and Shapiro [37]. This model, which interpolates between Wigner- 



Dyson statistics and Poisson statistics, is defined by the joint eigenvalue probability dis- 
tribution 

P(H)dH =dH J dUe ^+mH* + bTrUHUTH^ (1) 

where if is a Hermitian nxn matrix. The integral is over the unitary group with invariant 
measure denoted by dU. Critical statistics |36j is obtained in the thermodynamic limit 
with b scaling as b = h 2 n 2 at fixed h . In that case, the two-point correlation function 
decays exponentially at large distances and the number variance has an asymptotic linear 
behavior with slope less than one. In the thermodynamic limit, Wigner-Dyson statistics 
is obtained for a weaker n-dependence of b, and Poisson statistics is found for a stronger 
n-dependence of b. 

In this article we are interested in ensembles of nonhermitian random matrices. The 
study of random matrices with no restrictions imposed was initiated by the classical work 
of Ginibre [|lj. He found closed expressions for the two-point correlation function of the 
eigenvalues of a Gaussian ensemble of random matrices with complex entries. 

An ensemble that interpolates between the Ginibre ensemble and the Wigner-Dyson 



ensemble of Hermitian matrices was introduced in [30, 31 



P(C)dC ~ (fCe^ TrCtC+ ^) Tr(c2+(Ct)2) . (2) 

Here, C is an arbitrary nxn complex matrix with integration measure given by the 
product of the real and imaginary parts of the differentials of the matrix elements of C. 
For r = this model reduces to the Ginibre ensemble whereas for r = 1 (—1) it reduces 
to a Gaussian ensemble of ( ant i-) Hermitian matrices. The eigenvalues of this ensemble 
are scattered inside an ellipse with eccentricity given by 2-y/F/(l + r). 

The joint eigenvalue distribution can be obtained by using two alternative decompo- 
sitions 



C = UTU j and C — VAV~ l , (3) 

where U is a unitary matrix, V is a similarity transformation, T a upper triangular matrix 
and A a diagonal matrix. The diagonal matrix elements of T coincide with the complex 
eigenvalues A^k = z k- The invariant measure factorizes as (2[ 

dC ~ dU dTA({A kk })A({A* kk }) (4) 
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with the Vandermonde determinant defined by 



a ({^}) = no**- 2 *)- 

k<l 



(5) 



Since the Gaussian integral over the off-diagonal matrix elements of T factorizes it can 
be performed trivially. The integral over U is equal to the group volume. The joint 
probability distribution of the eigenvalues is thus given by 



P(A)rfA ~ dA |A(A)| 2 e~^^=° [N 



2 -f(4+K) 2 )l 



(6) 



This model has been analyzed in two domains: weak nonhermiticity and strong nonher- 
miticity. In the first case the thermodynamic limit is taken at fixed n(l — r), whereas in 
the case of strong nonhermiticity — 1 < r < 1 remains fixed for n — > oo. The two-point 
correlation function of this model was derived in |30] . 

In this article, we analyze a model that interpolates in between the models defined in 
eqs. (|l]) and (§). Our random matrix model is defined by 



P(C)dC ~ dC e^ c ' c -^^ 2+ ^ J dUe a * 



Tr UCU^C^ 



(7) 



where is C an arbitrary complex nxn matrix, and dll is the Haar measure of the unitary 
group U(n). In the special case of C being a normal matrix ([C, C'] = 0), a unitary 
transformation brings C to a diagonal form and the integral over U is the standard 



Itzykson-Zuber integral |47j given by 

dUe a 3 TrUCU^ 



det e 



A({z k })A({zl}Y 

where the Z( are the eigenvalues of C. One thus finds the joint eigenvalue distribution 



(8) 



P({z k }) ~ e -^"-[ ai|2l|2+ ^ (z ' + < 2) ]det 



(9) 



In the next paragraph we will show that this result is valid even if C is an arbitrary 
complex matrix that can be decomposed according to 

We start from the triangular decomposition C = UTU\ Since T is an upper-triangular 
matrix, the exponent in the integral over U in ([!]) is then given by 



(10) 



i<l 



After performing a trivial U(l) integration, the integral over U in (||) is over SU(n). The 
generating function for such integrals is given by 



/ dUe Tr( ~ Jut+J ' u) = Ffdet J, det jt, {Tr (J f J)H), 

JC/eSU(n) 
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where J is a complex nxn matrix and the functional form of the r.h.s., with k running over 
all positive integers, follows from the invariance of the group integral. In the expansion of 
the exponent ( |10|) all terms have the same number of factors U and U*. By differentiating 
fly]) with respect to J and J* at J = 0, we find that such terms can be only non vanishing 
if the sum of the indices of U is equal to the sum of the indices of U* (for the terms that 
enter in the expansion of the determinant the sum of the first indices is equal to to sum of 
the second indices). We thus find that in the expansion of (|H]) all terms with off-diagonal 
elements of T or T* vanish after integration. We conclude that the result (H) for the 
Itzykson-Zuber integral is also valid for an arbitrary complex matrix C with eigenvalues 

Zk- 

For convenience, the constants in the joint eigenvalue distribution of (|7|) will be pa- 
rameterized as 



A 



a 2 



1_ T 2> 

Ar Aa 2 



1 — r 2 r(l — a 2 ) ' 
r(l -a 2 )' 



«s = -r r^v ( 12 ) 



After a rescaling of the matrix elements of C by a factor 1/y/X the joint eigenvalue 
distribution of the model (0) reduces to 



P(A)dA ~ dAe ^H 1 - 21 11 "(i-tV 1 ar(i-«^ . • 'J det [e^=3»)^ . (13) 

We will analyze this model in two limits. The case when 1 — r remains finite in the 
thermodynamic limit will be referred to as strong nonhermiticity. In this class of models 
we will consider the limiting case of zero eccentricity 

a 

a^O, t -> with - = b fixed, (14) 

T 

which reduces to the Ginibre model in the limit in which the parameter b is taken to zero. 
On the other hand, the case of weak nonhermiticity |3(J is defined by the limit 



n — > oo, (1 — r)n = a 2 fixed. (15) 



Finally, let us mention that the wave functions of our model are distributed according 
to the invariant Haar measure of U(n). It could be that for diagonal U in (0) the wave 
functions show a multifractal behavior, but that this property is obscured by averaging 
over all U whereas eigenvalue correlations remain unaffected. 
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3 Strong nonhermiticity 



In this section we consider the case of strong nonhermiticity ([14]) . In order to rewrite 
the Itzykson-Zuber determinant in Eq. (p~3D in terms of an expectation value of two Slater 
determinants, we expand the exponential as 

oo jjm 



^ = E^"K) m - (16) 



m=0 

By a series of elementary manipulations we find 

oo oo tmi+...+m„ 

dete 6 ^ = E ••• E . ■ E (-l) a(7r) ^«(i)) mi • • • C"(4(n)) m " 

mi=0 m„=0 " ' 7reS„ 

£ 6 m 1+ ...+m n ^ (17) 

mi<m 2 <-<ra„ \J m r V m l- 

Including the other factors of the joint probability distribution we thus find 

P{z)dz~ E 6 mi+ - +m "det0 mj (^)det</> mi (4), (18) 

mi<ma<- <mn 

where the normalized wave functions are given by 

4> k (z) = -^r/^ 2 ' 2 (19) 

V 71" ft! 

satisfy the orthogonality relation 

#z<i>* k (z)<i> l (z) = 5 kl . (20) 

They are the single particle wave functions of the lowest Landau level of a particle with 
unit mass in a constant magnetic field perpendicular to the plane. The Hamiltonian of 
this system is given by (z = x + iy) 

H= l -{id x - y y + \{id y + xY. (21) 

and the corresponding Schrodinger equation reads 

HM?) = M?)- (22) 

If we write 

b = e~P (23) 

the joint probability distribution is equal to the diagonal element of the n-body den- 
sity matrix of the lowest Landau level fermions at temperature 1/(3, with an additional 
degeneracy-breaking Hamiltonian given by the absolute value of the angular momentum 

L = iyd x - ixd y , (24) 
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or equivalently of 



H = H + 2L = l -{id x + yf + l -{id y - xf. (25) 



The average spectral density p n (z), which can be interpreted as the one-particle density, 
is obtained by integrating the joint eigenvalue density over all coordinates except one. By 
using the orthogonality relations ( |Hf ) one easily finds 

1 n 

Pn{z) = — £ ^e'^-^^iz^iz*), (26) 

mi<ra 2 <-<m„ i=l 

or in an occupation number representation 

Pn(z) =y £ e~^ P pn P £ nfc ^(z)^(z*), (27) 

where the occupation number nk runs over {0, 1}. The partition function Z n is defined 
in the usual way 



ni+n 2 +--=n 



£ e ~PE p Pn r (28) 



Such sums can be easily evaluated in the grand canonical ensemble 

A n 



where we have introduced the prekernel 



E^§eV Z ), (29) 



{z x z? k 



H^)^-^ (30) 



" = £ n i + c-w (31) 



The fugacity ( is determined by the normalization of the one-particle density 

1 

For /3 <C 1 the sum can be converted into an integral resulting in 

C = e n/3 - 1. (32) 

Similarly, the two-point correlation function is obtained by integrating over all eigen- 
values except two. Again by going to the grand canonical ensemble one easily derives 
that the connected two-point correlation can be factorized in the result for the Ginibre 
ensemble and the prekernel (j30l) 

R2(z u z 2 ) = -\e-^- z ^\k(z 1 ,z 2 )\ 2 . (33) 



7 



For /3 C 1 but n8 >1 a partial resummation of the prekernel ( |3~0"D results in 

Ul ' " 2) " t Sfe! 4cosh 2 (/3(A;-n)/2)' (34) 

where x) = t fc_1 e~*(it is the incomplete T-function. For /3 — >■ it is justified to 
make the approximation 

11 8 

(35) 



1 + e P(k-n) i + e /3(fc+l-n) 4 cos h 2 (/5(fc - n)/2) ' 

In the remainder of this subsection we will evaluate the prekernel in several limiting 
situations. 

If the distance of z\ and z 2 (both inside the disk of eigenvalues) to the surface of the 
disk is much larger than 8, the numerator attains its maximum value when the Fermi- 
Dirac factor is close to unity. In that case the Fermi-Dirac distribution can be replaced 
by a sharp cutoff and the two-point correlation function is given by 

R 2 (zi,Z2) = -^e-^ 2 . (36) 

Inside the disk the average spectral density is The unfolded two-point spectral 

correlation function thus coincides with the Ginibre result. 

A more interesting situation arises in case both z\ and z 2 are close to the surface of 
the disk of eigenvalues. A nontrivial thermodynamic limit of the surface correlations is 
obtained for 

p ~ 4=, 



\z\z 2 \ ~ n, 

arg(zi^) ~ —=. (37) 



?? 



Using the asymptotic expansion for the incomplete T-function we find 



s P ^ Erfc((^* - k)/y/2k) 

[Zl,Z2) " V^h 4cosh 2 (/?(A;-n)/2) 



J3_ roc ^ Erfc((z 1 z*-n-t)/^2(n + t)) 

v^Fi-oo 4cosh 2 (/ft/2) ' ^ ' 

where Erfc(x) = f£° e~ t2 dt. We parameterize the vicinity of the surface of the domain of 
eigenvalues as 

Si ~\~ s* 

Zk = Vn + Sk, k — 1,2, and s = — - — -, (39) 
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where n>l and \s k \ <C \/n. Introducing the scaled temperature h by 



the prekernel simplifies for n —>■ oo to 



k(zi,z 2 ) = -j= dt -2- . 41 

^/7T J-oo 4 cosh t 



To the leading order in h, this expression can be simplified further 

k(zi,z 2 ) = -= f dye"* 3 f 



2 roc 2 f°° Q^Vlyht 

-= / dye~ y / dt- -5— 

V7T Jsv^ J-00 4 cosh t 



/"CO y p 

= 2v^/ dy^j—r. (42) 
For s ^> 1, the above integral dominated by the lower end point and is approximated by 

k(zi,z 2 ) ~ J- . 6 (43) 
V 2 sm(27rs/i) 

Accordingly, the spectral density near the edge to the leading order in h is given by 
p{z = \/n + s) = —k{z } z) = — -=■ / ay 



7T ' yfit Js sm(2iiyh) 

1 he~ 2s2 



^2~n sin(27rs/i) 



(44) 



At the zero temperature h — *■ 0, it reduces to the spectral density for the Ginibre ensemble 
close to the edge given by 0] p(s) = e~ 2s2 / (27r) 3//2 s. Likewise, the two-point function given 
by ( p3|) simplifies to 

1 /i 2 e - [( sl+s *) 2+ ( S2+s 2) 2 ]/ 2 

#2^1 = V™ + 51,^2 = V™ + «2) = " \ • I I 1 tU ^| 2 ' ( 45 ) 

27T I sm(7r(si + s 2 )h)y 

for |si + s|| ^> 1. As a consistency check, we find that the zero temperature limit for 
l/i - 2/2 > %k (with Sk = x k + iy k ) , 

1 e -2{x\+x 2 2 ) 

R 2 {z l: z 2 ) = -—- -. (46) 

2ir 6 (yi - y 2 ) 2 



is in agreement with the result in j|8| although different prefactors have appeared in the 
literature [49, 28 1. We mention that at zero temperature the asymptotic behavior of the 
prekernel can be obtained directly from its definition (|30| ) and agrees with (f46|). 
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On the other hand, in the high temperature limit the Fermi-Dirac distribution in (|30|) 
can be replaced by a Boltzmann distribution. The prekernel is thus given by 

k(z l ,z 2 ) = e-^Y, [ -^T L C^ k - (47) 
In this limit the fugacity is equal to ( = (3n, resulting in 

k(z u z 2 )=pn. (48) 
This requires us to define the scaled temperature by 

p = -k' (49) 

as opposed to the low-temperature case (flOp . The spectral density is thus given by 

p(z) = -1 (50) 
and the two-point correlation function has the exponential form 

R 2 (zi,z 2 ) = -^-e-^-^ 2 . (51) 

7T At 

Since the average spectral density decreases as 1/h, the unfolded eigenvalues become 
uncorrelated (Poisson statistics) in the high temperature limit. 

4 Weak nonhermiticity 

In the case of weak nonhermiticity, we start from the identity 



2, oo m 



e^**? = v/T^eMi^,^ ) £ " Hm ^) Hm {-L), (52) 

^ ml y/f ^ 

where H m (z) are the Hermite polynomials. Performing exactly the same manipulations 
as in (|17]) we obtain 



ft mH hm n 



mi<m2<-<m„ 

-.2 



x det det (53) 

ymj! V m i- 

The joint probability distribution (O) can thus be written as 

/nA miH hm " 

P(z) ~ 7T n (l - a 2 r/ 2 (l - r 2 )"/ 2 2 " det0 mi (^)det0 mfc (^), (54) 



mi<rri2<---<m 71 



T 
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where the wave functions defined by 

T k/2 



(55) 



satisfy the orthogonality relations [5C] 



d 2 ] z (j> k {z*)(j>i{z) = S ki . 



(56) 



The above wave functions ([55]) also span the set of the single particle wave functions in the 
lowest Landau level obeying the Schrodinger equation ( pT| - p2] ) which, in terms of properly 
rescaled coordinates, reads 



~{l-r 2 ){id x 



1-r 2 



;i-t 2 )k + t ^ 



If we write 



T 



(57) 



(5f 



the joint eigenvalue distribution may be interpreted as the diagonal element of the n-body 
density matrix of the lowest Landau level fermions at temperature 1/(3. The Schrodinger 



equation corresponding to fl25|) now reads 



1-r 2 



+ r(l-r)H 



x 



1-r 2 



+ 



r 



1-r 2 



{x + iyf 



<p m = (2m + 1) 



(59) 



Although, this relation is physically appealing we do not rely on it to obtain our results. 

Now we turn to the calculation of correlation functions. The p-particle correlation 
function is obtained by integrating P(zi, • • • , z n ) over z p+1 , ■ ■ ■ , z n . Using the orthogonality 
of the wave functions and expressing fl54|) as a single determinant one easily finds 



Pp (^1 j " " ' j Zp 



IV. 



(n — p)\ 
1 



d z p+ i ■ ■ ■ d z n P n (z) 



f E ,dft E' 

mi<m 2 <-<m„ l ' J '•••' P k=l 



-(3m k 



>rn k \Zi)<pm k \ z j 



(60) 



Here, the overall normalization constants Z n have been chosen such that the joint proba- 
bility integrates to unity. In an occupation number representation this correlator can be 
written as 



Z « n +ni+-=n M k=0 



k {Zi (p k {Z- 



(61) 
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where the occupation number rik runs over {0, 1}. Such sums are easily calculated in the 
grand canonical ensemble 



Rp\ z li ' i z p) 



j oo 



(62) 



n=0 



(63) 



where ( is the fugacity and Z is the grand canonical partition function given by 

oo 

^=n(i+ce- /3fc ). 

fc=0 

In the thermodynamic limit the correlators obtained by means of the grand canonical 
ensemble coincide with those from the canonical ensemble. The sum of the n& can now 
be performed easily. The result is given by 



Rp(zi, ■ ■ ■ ,Zp) = det K(zi,Zj), 
i,j=i,..., P 



with kernel defined by 



K 



(64) 



(65) 



The average spectral density, obtained by integrating over all eigenvalues except one, is 
thus given by 



p(z) = K(z,z) = ^ 1 + c _ le/3fc • 



The fugacity follows from the normalization integral and is given by 

1 



n 



(66) 



(67) 



fc=o 1 ^ •» e 

Similarly, the two-point correlation function is obtained by integrating over all eigenvalues 
except two. Subtracting p(z\)p(z-i) results is the connected two-point correlation function 
given by 

R2(zi,z 2 ) = -\K{ Zl ,z 2 )\ 2 . (68) 

As in the case of strong nonhermiticity, the kernel can be simplified by means of a 
partial resummation 

1 1 



K( Zi , Zj ) = Ys<t>k{Zi)<t>l{Zj) 

m=0 k=0 



1 _|_ £-l e /3m I + £-lg/3(m+l) 



m=0 

where the zero temperature kernel is defined by 

m 



1 _|_ ^-l e /3m X + ^-l e /3(m+l) 
1 1 



(69) 



(70) 



fe=0 
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4.1 Correlations in the bulk 



The bulk scaling limit of the zero temperature kernel fl70|) was analyzed in detail in |32[ 



We will recall their method for the sake of completeness. Using an integral representation 
of the Hermite polynomials, it can be rewritten as 

^0/ x 1 -^rM?M4\-¥A+4+z?+z?)]+^l+z?) 



x 



oo poo 



-oo J-oo ml 



1 

-P 2 



TT 2 T\/l — T 2 

roo roo 



( i ^ 3 j[|*i| B +|^|-i(a?+ai+*J 2 +a3 a )]+^(^+*5 2 ) 



x /" /"^^■c M2(1 ~ 1/T) ~ ,;2(1+1/T)+ ^ (zi ~ z * )/T+ ^ (2l+ ^ )/T r ^ m + 1,u2 - ^ 



i-ooJ-oo ml 
where r = u + v and s = u — v. The f-integral can be performed by a saddle-point 
approximation. To the leading order, the argument v in the incomplete T-function can 
be replaced by its saddle-point value given by 

For u 2 — v 2 ~ m and m — > oo, the incomplete T-function can be approximated by a step 
function 

1 x 2 

— -T(m + 1, u 2 — v 2 ) ~ 1 for u 2 < m + v 2 = m — — (73) 

ml (1 + r) 2 

and zero otherwise, depending on whether its integration domain contains the saddle 
point or not. We thus find the kernel 

\J m+v 2 



K m ( Zl ,Z 2 ) = / " due u\l-lH+iu( Zl -Z* 2 )/T 

tt 2 tV1 ~ ry 1 + 1/r J-V^+^ 

In the limit of weak nonhermiticity we magnify the bulk of the spectrum according 

Zi = xsjn + — — + i—, 
r nr y 2 

r 2 = 1--, (75) 
n 

where —2 < x < 2. For n — > oo this results in 



K°Jz u z 2 ) = J^e-^^i+i'^) / ^ )/B rf,» fi -^to)) (76) 



r K 2 a\[2 J-y/(m+v 2 )/n 
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For (3 — > the sum over m can be replaced by an integral. In this limit the kernel 
]) is given by 



K(z 1 ,z 2 ) 



i+x 2 /i" 4cosh 2 (/3nt/2) 



ndt 



n 



dt- 



2(3 



■s/l+t-x 2 /A 



7raV27r J-i+x 2 /4 4cosh \(3nt/2) Jo 

1 . 5 



due 2 cos u(irr + i(yi + y 2 )) 



xe -jj(y 2 i+y 2 2 )+^x( yi -y 2 ) 



n 



nay 2it 



where the combination 



dp 



1 _|_ e (p 2 -l+2 2 /4)//i 



e 2 e 



jp(7rr+i (i/i +3/2 ) ) „~ j*£ {v\ +V 2 ) + 1 31 (f i ~ V2 ) 



(77) 



(78) 



is kept fixed in the thermodynamic limit. Finally, we derive the small h limit of the kernel 
for x in the center of the spectrum (x ~ 0). The second integral in (|77| ) is rewritten by 
expressing the Gaussian term as, 



e 2 



: -~^+m(7rr+i(2/i+j/2)) 



1 



av 2vr «/-c 

After performing the integral over u we obtain 



dse 



(s-Tvr-i(y 1 +y 2 )) 



(79) 



K(z 1 ,z 2 ) 



n 



ir 2 a 2 



dse' 



-(s-irr-i(y 1 +y 2 )) 
^ 



-l cosh 7T7- 



+1& 



(80) 



The integral over t can be performed to leading order in h. In that case a/1 + 1 can be 
expanded to first order in t and the resulting integral over t, after extending its lower 
limit to — oo, is known analytically. We finally obtain 



K(z 1 ,z 2 ) 



nh 
2^? 



-e i 2 



dse' 



-(s-nr-i(y 1 +y 2 )) 



2a 



sin s 



sinh(irsh/2) 



(81) 



Sometimes it is useful to explicitly display the h = contribution to the kernel. From 
the second integral in Q77D at h = one can explicitly find the zero temperature result 
reported in |HJ. By subtracting and adding this term to ( 51 ) we find 



K(z l7 z 2 ) 



2n 1 

7TO \p2-K 



^l+«2 

=e a 2 



(an) 2 7lJl 

due 2 cos(w(7rr + + 2/2))) H 7= (82) 



sin s 



sin s 



sinh(7i7is/2) nhs/2^ 
where the first and third integrals cancel each other. 



e 2a 



2aV2vr 



±2-(s-(r7r+i(j/i+i/2))) 2 
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The spectral density at the center of the band is given by. 

2n 1 _ -'«- 



P{y) = K ( z = Wl Vn,z = iy/y/n) 



ira V27T 



pi a 2 t 2 

/ dte ~ cosh(2t?/) + 
Jo 



dt 



sint 



irh 
a\/27r 

sint \ _ t 2 2 
1 e 2^ cos(2yt/a ) 



s sinh(7r/it/2) 7rM/2 y 
where y\ = yi = y. The integral over lm(z) of the spectral density is given by 

/OO "j POO 

K(z,z)dlmz = —= p(y)dy 
-oo v 7T< J —oo 



n 



7T 



In Fig. 1, we show the normalized kernel K(z\, z?) defined by 

K(z 1 ,z 2 ) 



K(z 1 ,z 2 ) 



yjp(zi)p(z 2 ) 



(83) 



(84) 



55) 



for h = 0.15 and different values of the nonhermiticity parameter. We find that the 
spectral correlations weaken for increasing values of a and approach the result for the 
Ginibre ensemble for a w 2. Although not shown in the picture, it was verified numerically 
that the exact result ([77]) is almost indistinguishable from the small h result (^) for values 
of h up to h ~ 0.3, and significant differences are only found for values of h as large as 
The normalized critical kernel for Hermitian ensembles |37| is easily reproduced 
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from the ratio (]85f) starting from the expression ( |31~D and taking the limit a — > 0, 

7r/i sin(7rr) 

%^2H T . W2 L 86 

2 smh(7r"T/i/2) 



If we consider the a — ► limit of the kernel fl82|) or the spectral density fl83|), 5-functions 
of the imaginary part of the eigenvalues have to be taken into account carefully. For 
example, the a — ► limit of the spectral density (RBf) is given by 



71 

p{z) = -6(y). (87) 

7T 

Finally, let us mention that for a > 1 we recover the Ginibre's kernel for general 
complex matrices. 

4.2 Correlations at the edge 

Next we consider a microscopic scaling limit at the vicinity of either edge of the band of 
eigenvalues for z ~ ±2y/n, as an extension of edge correlation of the Hermitian Random 
Matrix ensembles. 

We shall need a more refined asymptotic formula for the incomplete T-function than 
Eq. (|73|). For x^m and m 3> 1, the incomplete T-function is dominated by the contribu- 
tion from the lower end point, so that pTI, |27J 



X 

r(m + l,x) = e" a 



x — m 



^ / m 
1 + 



(x — my 



(88) 



Accordingly, the kernel at zero temperature (|71|) reads 

K° m (z u z 2 ) ~ 1 p -j(T^)+il 2 +l^l-f(% 2 +4+4 2 +4 2 )]+^(% 2 +4 2 ) 



f°° r 00 dr ds 



x / / "" ^{-r 2 /2+irzi-s 2 /2-isz*)/r+{m+l) logrs-logm! (89) 

7-oo 7-oo rs — m 

For 21,22 ~ 2-y/n, m ~ n, and r ~ 1, the two saddle points of the r (s) integral merge 
at r = (s = —iypa). In order to obtain a nontrivial result, we magnify this region 
according to the scaling 

77i = n + n 1/3 t, (90) 

r 2 = 1--, 

n 

and change the integration variables as 

r = i\/n + n 1 ^p, s — —i>/n — n l ' & q. (91) 
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The subleading terms in ( jS8|) are of order 0(n 1 ' 6 ) in this scaling limit and can be ignored. 
To the leading order in n we obtain 

2^ 1/3 Xn^-JLM+vZ) f 00 r dpdqe4+ip(^)+4+i q ( X2 -t) 



V 71 a J-oo J-oo 



2tt 2tt — i(p + q) 

'l??!ljhi-V2)-$<rt+vZ) f dt > Ai(xi - t')Ai(x 2 - t'). (92) 
71 a J-oo 



where Ai(x) is the Airy function 

Ai(x) = r P 4+ ' PX = f - «»( £ + P*) ■ 03) 

./— rv-, V.TT In 7T \ .i / 



2tx Jo 7i 

The integral in Eq. (|92l) is called the Airy kernel K^{x\ ~ t, x 2 — t) (See Ref.0, §18), 
describing the edge correlations of the Gaussian Unitary Ensemble. By partial integrations 
one may express it in an alternative and more familiar form 

K Ai (x u x 2 ) = Ai(Xl)Ai/(X2) ~ Ai/(Xl)Ai(X2) . (94) 

Xi - x 2 

The scaling of m in (^) requires the introduction of a finite temperature parameter 
h by 

in contrast to the bulk scaling (|78|) . After replacing the sum over m by an integral over 
t, the low-temperature limit of the kernel fl69|) is given by 

K(z u z 2 ) = JlV—S y ^-^ + ^ r dtK M (x 1 -t,x 2 -t)^(—^- 7K 
\ Ti a J-n 2 / 3 dt \ 1 + e'/ ft 

'l^eiGtt-w)-^^) f°° it Ai(xi-t)Ai(x 2 -t) 

Due to the different orders of the level spacings in real and imaginary directions, the zero- 
temperature kernel is factorized, unlike the bulk kernel Eq. (52) of |30| or our Eq. (|76|). 
Namely, the dependence of on the order m is merely to dilate the eigenvalue support, 
which can be compensated by a change of the real part of the eigenvalue coordinate, 
x — > x—t. Accordingly, the effects of nonhermiticity and finite temperature are factorized. 
The former is reflected in the scaled kernel as a Gaussian blurring in the y-direction 
whereas, as the temperature h increases, the oscillation of the scaled spectral density along 
the x-direction is weakened toward the Poissonian limit. This is shown in Fig. 2 where we 
plot the spectral density in the Hermitian limit given by p(x) = J dtAi(x-t) 2 /(l + e t/h ). 
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5 Number variance 

The number variance in an arbitrary domain A of the complex plane is given by, 

S 2 (L) = L — [ d 2 Z! [ d 2 z 2 Y 2 (z u z 2 ) with L= [ d 2 zp(z), (97) 

J A J A J A 

where p(z) = K(z, z), Y 2 (zi, z 2 ) = \K(zi, z 2 )\ 2 and K(zi, z 2 ) is the spectral kernel defined 
in (|77|) . Apart from edge correlations we have found that in the strong nonhermiticity 
case the two-point correlations decay exponentially on a scale of one level spacing or less 
which results in an asymptotic linear dependence of the number variance on A with unit 
slope. Below we focus our analysis on the more interesting weak nonhermiticity limit. 

As will be seen in the figures below, the fluctuations of the eigenvalues increase with 
both increasing temperature h and increasing degree of week nonhermiticity a. The 
reasons for such behavior are the following: For larger values of h, the correlations of 
distant eigenvalues are suppressed resulting in stronger fluctuations and the slope of the 
asymptotically linear number variance increases with h. By increasing the degree of 
nonhermiticity, eigenvalues have more room to avoid each other along the imaginary axis. 
As a consequence, spectral fluctuations are stronger and deviations from Wigner statistics 
are observed. 

In the limit C 1 we calculate the number variance for the area A = [—L x /2, L x /2] x 
(— oo, oo). Because of the normalization integral (|8~4D we choose L x = Ln/y/n so that the 
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h = 0.1,a = 

h = 0.1,a= 1.5 

h = 0.1,a = 2.5 




Figure 3: The small h behavior of the number variance S 2 (L) versus L given in (|9~9"D for 
h = 0.1 and values of the nonhermiticity parameter as given in the legend of the figure. 



area A contains L eigenvalues on average. The dependence of the kernel on x is subleading 
in the thermodynamic limit. This allows us to rewrite the number variance as 



(98) 



z7r r lj r°° r°° 

E 2 (L) = L — / dr(nL/^/n-irr/y/n) / dy 1 dy 2 \K(z 1 , z 2 ) | 2 , 

n 6/z JO J-oo J-oo 

where the prefactor includes a contribution from the Jacobian of the transformation (|75|). 
The integrals over yi and y 2 are easily performed in terms of the variables u = y\ + y 2 and 
v = y\ — y 2 . The final result for the small h limit of the number variance is thus given by 



TX 2 h 2 



S 2 (L) = L-2 f dr{L-r 
Jo 

sin 2 (t) sin 2 (t) 



sin 7rr 



7r 2 r 2 



+ 



4 



dt 



smh 2 (7iht/2) {nht/2)' 



_ 1 a-nr) 2 



(99) 



We observe that in this limit the finite temperature effects decouple from the weak 
nonhermiticity corrections. For L ^> 1/h and a L it can be shown from fl9"9] ) that the 
number variance is given by 



E 2 (L) 



7T 



3/2 



^ + ^L + 0(1/L), 



7T Z 



(100) 



where 7 is the Euler constant. The term linear in a can be calculated in the h — * limit 
and was obtained in |J0|, whereas the term linear in h can be calculated for a — > and 
was derived in p7| . In Fig. 3, we show the small h limit of the number variance (|99|) 
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6 - 



— ' 1 ' 

— h = 0.1, a = 20 
- h = 0.1,a = 25 
-■ h = 0.1,a = 30 




Figure 4: The number variance 
parameter a. 



is computed for large values of the nonhermiticity 



for h = 0.1 and different values of the nonhermiticity parameter. We observe that the 
asymptotic linear behavior given by (|100| ) is already reached well below the expected scale 
of 1/h. We remark that for values of h as large as 0.3 the small h result (|97D is still very 
close to the exact result obtained with the kernel (|77|). 

The small h result for the number variance (^9]) is also valid for large values of the 
nonhermiticity parameter. Plots of (|99| ) for a ^> 1 are shown in Fig. 4. We find that the 
asymptotic result for the slope is still approximately given by h/2 and depends only weakly 
on a. For L < a we find that T?(L) — > L which is the result for strong nonhermiticity. 
This crossover behavior was first found in the limit h — > p0| . 

The imaginary part of the eigenvalues is of order a. This is shown in Fig. 5 where 
we plot the p(y)/p(0) (with p(y) given in eq. (|83f) ) versus y. Since the imaginary part of 
the eigenvalues is of the same order as the spacing of the real part of the eigenvalues, the 
number variance computed for a rectangle < Imz < Ay C a is expected to be given by 
S 2 (L) — ► L where L is the total number of eigenvalues in the rectangle. This is shown in 
Fig. 5 where we plot the number variance obtained from 



using the kernel fl82|). 



6 Conclusions 



In this article we have introduced a two parameter ensemble of complex random matrices 
with no hermiticity conditions imposed. This ensemble interpolates between the Gaussian 
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Figure 5: The renormalized spectral density p r {y) = p(y)/p(0) (with p{y) defined in (|83|) 
in the center of the band is shown for different values of the nonhermiticity parameter. 



Unitary Ensemble, the Ginibre ensemble and the Poisson ensemble. Using methods from 
statistical mechanics and properties of orthogonal polynomials, we have analyzed this 
ensemble in two different limits: weak nonhermiticity and strong nonhermiticity. 

We have shown that the joint eigenvalue distribution of our random matrix model 
coincides with the diagonal element of the density matrix of a two dimensional gas of 
spinless fermions in the lowest Landau level at finite temperature. The two parameters 
of our model have been interpreted in terms of a shape parameter of the two dimensional 
domain of eigenvalues (or particles) and a temperature. 

In the strong nonhermiticity limit, in the bulk of the spectrum, the correlations of 
the eigenvalues are given by Ginibre statistics and decrease exponentially on the scale 
of the average level spacing. The situation is different near the surface of the spectrum 
where, at zero temperature, the correlations decrease as an inverse square law in the 
direction of the surface. At finite temperature this power-law behavior changes into an 
exponential behavior. At very high temperatures the surface and the bulk are no longer 
distinguishable. In that case the two-point correlation function of the unfolded eigenvalues 
still decays exponentially but with an exponent that is proportional to the temperature. 
In this way the Poisson limit is recovered at high temperatures. 

In the weak nonhermiticity limit there is no clear distinction between bulk and surface 
and the temperature affects the correlation functions of the eigenvalues. In the low tem- 
perature limit we have obtained a closed analytical expression for the two-point correlation 
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L 

Figure 6: The number variance given by the general formula fl97|). The domain of inte- 
gration is a rectangle in the complex plane containing L eigenvalues and width given by 
< \mz < Ay. The nonhermiticity parameter is equal to a = 0.4 and the value of h is 
equal h — 0.1 for all curves. The number variance is almost Poisson for Ay < a. 

function which reproduces critical statistics. We have found that, although level repulsion 
is still present, the number variance is asymptotically linear with a slope depending on the 
temperature parameter but not on the nonhermiticity parameter. A remarkable feature 
is that temperature and weak nonhermiticity effects decouple in this region. Thus critical 
statistics is not modified by a weak nonhermitian perturbation. 

Finally, let us explain a physical prediction of the present model. Since for critical 
statistics the slope of the number variance is related to the multifractal dimension of 
the wave function and, in our model, the slope does not depend on the nonhermiticity 
parameter, we predict that the multifractal dimension of a physical system does not 
depend on the nonhermiticity parameter either. We thus predict the same multifractal 
dimensions for open and dissipative systems. A simple model for which this prediction 
may be tested is a three dimensional disordered system at the critical density of impurities 
and with several leads attached to it. We thus expect that in the weak nonhermiticity 
domain the leads do not affect the multifractal dimension of the wavef unctions. 
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